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Abstract 

Sodium impurities are diffused electrically to the oxide-semiconductor interface 
of a silicon MOSFET to create an impurity band. At low temperature and at low 
electron density, the band is split into an upper and a lower sections under the 
influence of Coulomb interactions. We used magnetoconductivity measurements to 
provide evidence for the existence of Hubbard bands and determine the nature of 
the states in each band. 
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1 Introduction 



In the early 70's a substantial effort was put into the study of the effect of 
impurities in silicon MOSFETs. The main reason for this was to understand 
the origin of the instabilities in field effect transistors and the umpredictable 
changes in the threshold voltages. These studies contributed to improve the 
performance of a new generation of silicon-based chips for computer appli- 
cations and their reliability. In fact, the presence of impurities or defects in 
the silicon oxide leads to an inhomogeneity in the electron distribution at the 
Si-Si02 interface and traps electrons. The most common impurities are alkali 
metals like lithium, potasium or sodium. These trap states give rise to the 
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formation of an impurity band below the conduction band. This effect was 
first observed by Fang, Hartstein and Pepper [TH21I3] for high impurity concen- 
trations (> 1 x 10 12 cm" 2 ). A consequence of the existence of such an impurity 
band is that the onset voltage for conduction in oxide-doped MOSFETs is also 
shifted and the electron mobility substantially decreased. The impurity band 
produced by such a doping can be described by tight-binding models for high 
impurity density [I] whereas Mott-Hubbard models [5] need to be used for the 
low concentrations for which Coulomb interaction plays a defining role. Given 
single valency atoms like sodium, it would be expected that the electronic 
states in the band are made of single trapped electrons. Under the influence 
of electron-electron interactions, it has been suggested however that a stable 
state with two bound electrons would exist and would be characterized by a 
long lifetime (D_ state) [6]T7] . The band formed by the D_ states is referred 
as the upper Hubbard band (UHB) and the one formed with neutral states 
as the lower Hubbard band (LHB). The question of the existence of Hub- 
bard bands in Silicon MOSFETs has been put forward by Mott nearly thirty 
years ago to explain the magnetoconductivity of short disordered MOSFETs 
but the bands were never directly observed in experiments j8]. The study of 
Hubbard bands in semiconductors has regained some attention since the end 
of the 90's with the development of quantum information and quantum com- 
putation [9]. Following the difficulty in reading out directly the value of the 
spin in architectures developed in agreement with the Kane model, a spin to 
charge conversion was proposed. In this approach a stable D_ state is used 
to read out the result of the quantum operations. Some optical studies have 
previously been carried out in Mott-Hubbard insulators like Sr 2 Cu03 [TU] or 
boron-doped diamond [11] but no direct conductivity measurement has been 
performed so far. In the present paper, we have used magnetoconductivity 
measurement to provide direct evidence of the presence of Hubbard bands in 
a sodium-doped silicon MOSFET. We also describe the electronic configura- 
tion of the obtained D_ state by distinguishing between Is 2 , ls2s and ls2p 
states. 



2 Experiment 



All measurements were performed on silicon MOSFETs. Such devices have 
been widely used because of the ability to continuously vary the electron den- 
sity and the Fermi energy by use of a metal gate. The geometry of the device 
was chosen to be circular (Corbino geometry) to avoid leakage current paths 
around the source and drain contacts. The devices were fabricated using a high 
resistivity (10 4 f2.cm) (100) p-silicon wafer to minimize, as much as possible, 
the scattering with boron acceptor impurities, especially close to the silicon- 
oxide interface. A 35 nm thick gate thermal oxide was grown at 950 °C in a dry, 
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Fig. 1. a) Cross-section view of a corbino MOSFET used in the experiment when the 
sodium ions are close to the Si-Si02 interface, b) Schematic diagram of the density 
of states (DOS) for the present device, with a low energy (LHB) and a high energy 
(UHB) impurity band separated by a gap to the conduction band (CB). 

chlorine-free oxygen atmosphere. The sidewalls of the oxide were protected by 
thick insulating LOCOS (Local Oxidation of Silicon). The effective gate length 
of the Corbino MOSFETs was 1 fj,m and the diameter of the interior contact 
was llOyum. Contacts were realized by implanting phosphorous at high dose 
and sputtering aluminium. The contact resistivity was measured to be 3.5 and 
2.3f2.cm _1 at nitrogen and helium temperatures, respectively, and the sheet 
resistance was 6.3 and 5.9 fi.D" 1 for the same temperatures. Sodium ions were 
introduced onto the oxide surface by immersing the device in a lCr 7 N solution 
(~ 6.4 x 10 11 cm -2 ) of high-purity sodium chloride in deionized water. The 
surface of the chip was then dried with nitrogen gas and an aluminium gate 
subsequently evaporated. The application of a positive gate voltage (+4V at 
65°C for lOmins) causes the sodium ions to drift towards the Si-Si02 interface 
while the application of —4 V DC in the same conditions removes the ions from 
the interface. The ions are frozen at their position once the device temperature 
becomes lower than approximately 150 K (Fig. 1). Standard low- noise lock-in 
techniques with an amplifier gain of 10 8 V/A were used to measure the source 
to drain conductivity. An AC excitation of amplitude Vac = 15 fiV and a fre- 
quency of 11 Hz were chosen. The DC offset of the amplifier was suppressed 
using an appropriate blocking capacitor. The gate voltage was controlled by 
a high resolution digital to analog converter and the temperature measured 
by a calibrated germanium thermometer. The magnetic field was produced 
by an Oxford 12 T superconducting magnet and applied perpendicular to the 
Si-Si02 interface. 

Several devices were processed identically and gave results that lead to iden- 
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tical conclusions although we noticed some variation in the relative position 
and width of the impurity bands, as well as in the conductivity values. From 
the difference in the threshold voltage at 77 K, we obtained an effective ion 
density of ~ 3.7 x fO 11 cm -2 at the interface indicating that only 60% of the 
ions drifted to the interface. We also fabricated a number of control devices 
that were not exposed to sodium contamination and were used for comparison. 
The following results are presented for a specific device that was chosen for 
its high reproducibility in time as well as for its high signal to noise ratio. 



3 Results and discussion 

Fig. 2a represents the source-drain conductivity obtained at different values 
of the magnetic field. The dependence of conductivity on temperature for 
the same device in the hopping regime showed the presence of two groups of 
peaks clustered around V g = — 2 and —0.5 V. These were attributed to a split 
impurity band due to the presence of sodium impurities at the Si-SiC>2 interface 
[T2j . Arguments in favour of Hubbard bands were provided by studying the 
variation of conductivity at higher temperature, where activation mechanisms 
determine the behaviour [T3] . 



3.1 Negative magnetoconductivity 

The magnetoconductivity was found to be negative for the whole range of 
gate voltages studied. This is expected and was already observed in localized 
systems [141151116] (Fig. 2b). The variation with magnetic field is well described 
by In a ~ —aB 2 even at fields up to 5 T where a is gate voltage dependent. 
This behaviour is often attributed to an orbital compression of the donor 
wavefunctions due to the magnetic field. This results in a reduction of the 
wavefunction overlap between neighbouring hopping sites and an increase of 
the localization of electrons [17]. The wavefunction then acquires a phase factor 
(p ~ (^r 3 ) 1 ^ 2 jl\ where Ib is the magnetic length given by Ib = (h/eB) 1 ^ 2 , r is 
the hopping length and £ is the localization length. This assumes a gaussian 
distribution for the flux <fi. For transport processes involving only one hop 
within an area (£r 3 ) 1//2 as in the case of nearest-neighbour hopping or resonant 
tunneling, In a is found to be quadratic in magnetic field with a given by [IB] 




(1) 
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Fig. 2. a) Conductivity versus gate voltage at B = 0T, 2.5 T and 5T (from top to 
bottom) at 309 mK. For clarity curves are separated by a constant offset of 10 _3 e 2 //i. 
b) Magnetoconductivity for V g = -0.24 (o), -0.56 (•), -0.97V (*) and -2.08 V 

(<)■ 

Nt is the density of active traps and B c is a constant describing the number 
of bonds in a circular percolation problem, being tt in tunnelling problems and 
respectively 4.5 and 2.7 in the case of nearest-neighbour hopping in a zero or 
finite width impurity band. If, however, the transport is described by several 
hops, a minimization procedure, as first used by Mott, [19] needs to be carried 
out in order to find the optimum hopping length and hopping energy. This 
makes a temperature dependent [2D] . 



a 



1 eT 1 /T \ 3 (^) 



24 h 2 (p + 3) V T 



where p is in case of the Mott variable range hopping regime [21] or 1 in 
case electron correlations are important (Efros and Shklovskii regime [22] ) 



3.2 Low temperature conductivity and localization length 



The proper analysis of o~(B) is thus constrained to the correct description of 
the transport mechanism at low temperature. Below 1 K, the conductivity de- 
creases rapidly with temperature as Act = a (T) — a (0) ~ T 2,7 (Fig. 3a). This 
variation is not affected by the source-drain voltage and cVac < k^T in all 
measurements (Fig. 3b). Consequently, we do not expect that this behaviour 
results from electron heating in the contacts, but from a contribution of the 
acoustic phonons at low temperature. In our previous study, the same data 
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Fig. 3. a) Variation of the conductivity in temperature for = —0.24 (o), —0.56 
(•),— 0.97 V (*) and —2.08 V (<). b) Variation of the temperature exponent with 
Vac a t the same gate voltages, c) Variation of the localization length in terms of V g 
at 309 mK deduced from Eq. 2. 

was analyzed in temperature and we showed that the transport was charac- 
terized by hopping conduction with an exponent p ~ 0.4 above IK [12]. In 
the framework of 2D resistance networks [23] the conductivity is a ~ v/T, 
whereas the scattering time is r ~ ur 2 /T where v is the hopping rate and r 
the hopping length. Our results showed the best fits were obtained with an ex- 
ponential prefactor ~ j 1-0 ' 8 . This implies r varies as T~ L6 . This value is close 
to the T~ 3 / 2 law expected in non-polar semiconductors for acoustic phonon 
scattering [23] and as observed in silicon [25]. This excellent agreement with 
theory implies that acoustic phonon-assisted hopping is the main mechanism 
for transport above 1 K and that this transport is in fact still active and pre- 
dominant at 300 mK. Due to the short hopping distance, typically r ~ 3£, 
an electron is scattered by only a few impurities between the initial and the 
final sites. Consequently, Eq. 2 should be used with p = 0.4 to estimate the 
value of the localization length £ at low temperature, giving the values of To 
from our previous analysis [12]. Results are shown on Fig. 3c. The variation 
of £ (V^) is consistent with our previous observations, i.e. an increase in the 
localization at the band positions, but the values are now nearly half of those 
obtained from T Q (£) in the absence of magnetic field [12]. Such a halving of 
£ has already been observed experimentally in insulating devices with strong 
spin-orbit coupling [26J. It is plausible that such an effect takes place in our 
device because of the presence of strong potential fluctuations at the Si-Si02 
interface. We also observed that the conductivity of the lower band is sup- 
pressed at a field of about 3.5 T whereas the upper band is still conducting 
even at the higher field. This suggests that the lower band is energically deeper 
and is already strongly localized even in the absence of a magnetic field. Ef- 
fectively, the upper band is supposed to be closer to the extended states in 
energy than the lower band and it is likely to be more strongly influenced 
by the enhancement of the localization due to the shrinkage of the electron 
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wave function caused by the magnetic field. A point of interest is the absence 
of significant differences in the values of the localization length in the upper 
and lower bands. One may expect indeed the lower band to be more deeply 
localized and have a lower localization length. In fact, one needs to consider 
that, in our device and because of disorder, the lower Hubbard band may be 
formed by a majority of singly occupied states and a minority of empty or 
doubly-occupied states. In order to get transport through the lower Hubbard 
band, most of the electrons will hop to a singly-occupied state and form a D_ 
state. On the other hand, in the upper Hubbard band, most of the sites will 
be doubly-occupied and the electrons will hop directly to the conduction band 
edge because of the presence of a large conduction band tail and the presence 
of states in the gap between the upper Hubbard band and the conduction 
band. As a consequence, we expect to find the same localization length for the 
upper and lower band in transport measurement in the hopping regime. 



3. 3 Magnetoconductivity fluctuations 



We observed also a non-monotonic decrease in the conductivity with mag- 
netic field. This takes the form of fluctuations with amplitude dependent on 
gate voltage. These fluctuations are present in many devices except for those 
without sodium ions. The fluctuations appear as a series of positive and nega- 
tive magnetoconductivities. Their position in magnetic field also changes with 
time, but in general is periodic in B. Such fluctuations have already been ob- 
served by Nguyen [27] and were attributed to random interference between dif- 
ferent hopping paths. To precede with the analysis, the magnetic field value B m 
of the fluctuation extrema were extracted from the curve Aa(B) = o~(B) — <t 
exp(— aB 2 ) and plotted as a function of the extrema index (Fig. 4). We found 
a linear dependence with a slope B^ for all gate voltages indicating the fluctu- 
ations are periodic in magnetic field. To minimize the uncertainty due to the 
limited number of points, the experiment was repeated four times with the 
same configuration. We found that ^ = (h/eB^) 1 ^ 2 > £ but ^ ~ (r^) 1 ^ 2 . We 
can deduce the modulation corresponds to an integer number of elementary 
flux quanta through a surface of area 7rr£ and that B m is given by : 



This situation is similar to that giving rise to Aharonov-Bohm oscillations in 
quantum dot systems [28] with (r£ ) playing the role of the dot radius. In the 
hopping regime, the interference area has a cigar-shape with radius £ in one 
direction and the hopping length in the hop direction. Quantum interferences 
between forward and backscattering paths induced by the localizing impurity 
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Fig. 4. Measure of the critical field of the fluctuations and the corresponding values 
£ for Vg = -0.43 V, -0.56 V-1.50 V and -2.08 V at 309 mK. 

potential due to the magnetic field may be responsible for such fluctuations. 



4 Peak positions in magnetic field 



Finally, we discuss the variation in the position of the peaks in magnetic 
field (Fig. 5). We observed that the group of peaks clustered in the upper 
band moves towards more negative gate voltages with approximately a slope 
of ~ — 4mV/T whereas the one in the lower band move with a slope of ~ 
+4mV/T. Such movement in energy levels are expected in magnetic field, 
and the difference in the sign can be explained by including the effect of the 
orbital momentum. In fact, this behaviour cannot be interpreted as being due 
to spin effects as no spin splitting was observed. Also, the gap energy between 
the upper and lower bands estimated previously [13] is much larger than the 
spin splitting at 5 T. For simplicity, we suppose the electrons are localized at 
the impurity sites with an approximate parabolic potential of width uo. The 
variation AE of an energy level is then described by 



AE = (2n+ \l\ + l)h 



luJ c , 2 



-IfVjJr 

2 c 



(4) 



where n is the index for the Landau levels, I is the value of the angular mo- 
mentum and u c = eB/m* with m* = 0.19m e _ the transverse effective mass 
of the electrons in < 100 > silicon. 
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Fig. 5. Conductivity versus gate voltage and magnetic field at 309 mK. 

Consequently, n = and I = are associated with the lower group of peaks 
(non-degenerate) and n = 1 and / = 1 to the upper group (degeneracy of 2). 
The value I = 1 also implies the presence of two electrons in the configuration 
ls2p. This behaviour supports the idea of the formation of Hubbard bands [31] 
and a stable Na~ ion at the Si-Si02 interface as well as the presence of a strong 
spin-orbit coupling. This shows as well the importance of Coulomb interactions 
in this system. When increasing the magnetic field, the wave-function of the 
lowest state shrinks and the interaction between the two degenerate electrons 
grows until the Zeeman energy is sufficient to spin-flip the second electron 
and promote it to the I = 1 state. This transition is predominantly caused by 
Coulomb interactions since no spin splitting was observed and was still quite 
small even at 5 T. Finally, the ls2p configuration holds as long as AE < for 
the upper band, which is verified up to 5T. This implies that hu>o > 2.7 meV 
and provides a lower bound for the mean impurity potential depth. 



5 Conclusion 



We have observed a negative magneto-conductivity up to 5 T. This result con- 
firms previous observations on the same device that the system is strongly 
localized with a localization length of the order of the mean distance be- 
tween impurities. The fluctuations observed in the magnetoconductivity have 
been attributed to quantum interference due to the presence of impurity cen- 
tres close to the interface. Finally, the shift in the position of the bands are 
compatible with the presence of upper and lower Hubbard bands and a pos- 
sible transition from a singlet to a triplet state in the upper band induced by 
electron-electron interaction. 
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